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1 . 



Introduction 



The central issue in approximating a partial differential equation by 
a finite difference equation is the degree to which the difference equation 
solution agrees with the solution to the differential equation. This 
agreement can be considered in both its quantitative (e.g., relative 
error), or qualitative aspects, e*g., behavior of transients, propagation 
of fronts, etc. 

In this paper we investigate some of these aspects for a prototype 
dispersive wave (hyperbolic) equation. This type of problem is 
especially important in numerical weather prediction, since for ”fast” 
time scales, the meteorological primitive equations are dispersive. We 
consider both the differential and quasi-^discrete case, where second order 
centered spatial differences are used. In both cases solutions are 
obtained by Fourier transform methods, with steady state solutions 
extracted directly from the transforms and inverted in closed form, and 
the asymptotic behavior of the transients determined by the method of 
stationary phase. 

Our analysis shows that discretization of spatial derivatives has 
two major effects: 

(1) The initial conditions of the differential problem contribute 
to the steady state solution in a manner that decays exponentially with 
distance. This qualitative effect is retained in the quasi-discrete 
formulation, but the rate of decay is decreased. 

(2) The discretization introduces additional transients beyond 
those encountered in the differential case. These new transients are at 



2 



least of the same magnitude of decay as the differential transients, and, 



in some instances, they are more persistent. 
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2 . 



The Differential One'-Dimensional Adjustment Process 



One of the simplest models of dispersive waves is the linearized one- 
dimensional shallow water equations with no mean flow, in an infinite 
region: 



3u 

at 



fv + g 






|2 + fu . 0 



Ih + H - n 



( 1 ) 

( 2 ) 

(3) 



where u is the perturbation velocity in the x direction, v is the 
perturbation velocity normal to the x direction, H and h are the mean 
and perturbed heights of the free surface, respectively, and g > 0 and 
f > 0 are gravitational and Coriolis parameters, respectively. This model 
is especially important in the study of the meteorological problem of 
geostrophic adjustment, and has been studied in some detail by Rossby [1], 
Cahn [2], Blumen [3] and Winninghoff [4]. In their papers, the model has 
been studied by eliminating between the equations to arrive at: 



2 2 
3 u . .2 „ 3 u . 

r + f u - gH r = 0 , 



( 4 ) 



at 



3x 



then solving (4) by a Fourier Transform approach. After solving (4), 
solutions for h and v are obtained by substitution into (2) and (3), 
although closed-form solutions are not produced in some of the papers. 
Note the dispersive character of (4) is clearly seen by assuming a wave 
solution: 



4 



u(x,t) = A e 



i(kx-vt) 



which leads immediately to: 



+ k^gH = f^(l+x\^), X = /gH/f . 



( 5 ) 



We now derive an alternative means of solving (l)-(3) which is 
superior in that it does not require elimination, it produces u, v and h 
without back substitution, it yields interesting insights into the 
transient and steady state behavior of the solutions in the differential 
case, and it has an extension in the quasi-discrete case that is quite 
illuminating. 

We denote Fourier Transforms by a wavy overhead bar, e.g.: 



oo 




-ikx 



dx , 



( 6 ) 



etc. Then (l)-(3) reduce directly to: 




(7) 




( 8 ) 




(9) 



together with initial conditions; 



OO 




u(x, 0)e 



-ikx 



dx 



( 10 ) 



—CO 
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etc. Since (6)- (9) is a coupled set of constant coefficient ordinary 
differential equations, it can be solved by the usual process of finding 
the eigenvalues and eigenvectors of the coefficient matrix. This leads 



straightforwardly to: 



2 ivt 3 -ivt 
u == — e - — e 
V V 



, ±f ivt if -ivt . 

V = ikg ^2^ 2 ^3^ * ^ 

V V 



I:- . kH ivt kH -ivt 

h “ f “i - -2 “2e - ~ “3® 

V V 



( 11 ) 



where v is given by (5), and the are picked to satisfy the initial 

ivt ^ivt 

conditions. Observe that the e and e terms both represent the 

transforms of transients in the time domain. Solving the initial 



conditions for the , collecting like terms in (11), and simplifying 



yields: 



Or 

^ f V ikg h 

^ ^ o . o 

u(k,t) = u cos vt H sin vt 

o V V 



sin vt , 



f . . . fk^gH , . ikgf 

v(k,t) = — u sin vt + < — ° — I ;r cos vt > v H ^ 

VO § ^ ^ 2 

f V V I V 



|l - cos vt I , 



>( 12 ) 



2^.- V ikH 'v ikHf 

h(k,t )= — — u^sin vt — 



{i - + K 



. % 

^ cos vt>h 

2 i O 



Observe that (12) immediately yields by inspection the transform of 



the steady state (**balanced”) solution: 



us(k) 

Vs(k) 



?Js(k) 



0 



2o 2o o 2lfo 

V V V 



V 



} 



IkHf ‘i; 

— r- V + — r h 
2 o 2 o 

V V 







>v 



>(13) 



These can be immediately inverted by the convolution theorem, and the obser- 
vation that 




2X 



2X ^2 



2 2 
1+k X 



to yield: 



Ug(x) = 0 



Vg(x) = v(x,0) + 



hg (x) = h(x,0) 



00 

f 



-|x-s|/X . . 

e ^ (s>0) - v(s,0) ( ds 



>(14) 



2X^f 



oo 

I 



- I X-S I /X 



K ik 

f dx 



(s,0) - v(s,0) > ds 



Solving equations (l)-(3) by this approach allows us to easily 
observe several elementary results: 

(1) Since u (k,0) does not contribute to any of the steady state 
solutions, any "imbalance" in u(x,0) does not affect the final balanced 
state. 

(2) The term in (14) 



_ ^ ik 

f 3x 



(s,0) + v(s,0)/ 
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1 



the initial value of the quantity usually referred to as the ageostrophic 
wind, can be considered as the measure of "imbalance" in the initial state 
that is not dissipated at the final state. 

(3) At steady state, the effects of this initial imbalance contribute 
most strongly in the immediate neighborhood of the imbalance and die off 
exponentially in space away from it. 

Observe also that the transient part of (12) can be written: 



u^(k,t) 



% 

u cos V 
o 



t + — < V - fi ) sin vt 
V / o f o 






- — u sin vt + <1 V - } cos vt 

VO ^2)0 f o 



h^(k,t) = 



ikH 'V . ikHf ) 'v xkg y ( 

u sin vt z- < V ^ h } cos vt 

VO 2 ) o f o 






(15) 



Explicit inverses to these transforms are expressible in terms of 
convolutions involving J 

I 

the as}mptotlc behavior of the transients is possible by using the method 
of stationary phase. Let 






xW - 



). 



However, a simplified view of 



d(x,t) = v(x,t) - 1“ (x,t) . 



Then 



u^(x,t) 



f 

2tt J 



u(k,0)cos vt e^^^dx 






0)sin vt 



ikXj, 

e dk 



(16) 



Thus, using the method of stationary phase (Appendix 1), we can show, for 
fixed X as t : 



8 



r 






(x,t) 



2 2 
Af t 



2,(x2f2,2.,2)3/2 






il , , 0 



I /,2.2 2 2 

WX f t -X 



cos 






1 , 2^2 2 2 

\IX f t -X 

Xft 



JxW-x^ 



, 0 



. , 1 /, 2 . 2^2 2 . , 
Sin / y X f t -X 



( 17 ) 



where and ipy^ are slowly varying phases. When |x| < < Xft , this 

is more conveniently approximated as 



u^(x,t) 





2 


1 / \ 


1 




hi ^ 0 1 


2 




1 1 1 2 ’ ® I 


2'it X ft 


1 





cos (f t+<j)j^) 



+ 



(* ■ ”) I' 



2 , 0 I jsin(ft+i|;j^) 



. 



( 18 ) 



Thus the decay to steady state of u at a fixed point x is governed 

by two factors: 

- 1/2 

(a) A t decay, which can be interpreted as the effect due to 

the dispersive nature of the process, and 

(b) An additional possible decay, depending on the degree to which 
the initial imbalances were distributed in the long wave numbers, since 



« . / X 

iim u I — r — , 
t^ • \x^ft 



0^ I = iim u(s,0)e ^ 

oo 

■/ 



-i(x/X ft)s 



ds 



u(s,0)ds , etc. 
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Similar analysis of v^(x,t) and h^(x,t) yield 



v^(x,t) 



2tt X^ft 






, u , ,sin(ft+(fij^) 



+ 3 ' * 



(* 



, 0 



cos(ft+i|;^) y , 



( 19 ) 



and 



h^(x,t) 'V 



-I 

2 



2tt X^ft 



* te • •)! 



sin(f t+(j)j^) 



* K* ■ ")l 



COS (f 



( 20 ) 



Note the one difference in h^(x,t). Due to the presence of the (ik) 

term in h^(k,t), the decay due to dispersion of h^(x,t) proceeds as 

-3/2 -1/2 

t , rather than t 

In summary then, we have shown here that the differential formulation 
of the dispersive wave model (l)-(3) is most easily solved by Fourier 
transforming to a system of ordinary differential equations. This solution 
shows the model always tends to a steady state, whose difference from the 
initial state is determined solely by the initial ageostrophic wind field. 
The contributions from regions where the ageostrophic field is initally 
non-zero die out in the steady state exponentially with distance. Lastly, 

both the initial u and ageostrophic wind fields contribute to transients 

-1/2 -3/2 

that die out in time as t for the velocities, and t for the 

free surface height. 
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3. 



Second Order Centered Finite Differences - The Quasi-Discrete Adjust- 



ment 



In this section we present an analytic treatment of the most common 
difference scheme used for (l)-(3), and investigate the resulting effect 
on the steady-state and transient behavior discussed above. 

Consider the continuous, quasi-discrete, second order centered-leap-* 
frog formulation of (l)-(3) 



= f v(x,t) - -^[h(x+Ax,t) - h(x-Ax,t)] 
f =-f u(x,t) 

H = - 2 ^[u(x+Ax,t) - u(x-Ax, t)] 

If we Fourier Transform this system, noting 



y. (21) 



J 



00 

1 u(x+Ax,t)e dx 
— 00 



ikAx . 

e u(x,t) 



y 



we have: 



^ = f 

dt 



V 



_ ISl 
Ax 



sin kAx h 



dv 

dt 



Or 

- f U 



dh 

dt 



iH 

- — sin kAx u 
Ax 






y .( 22 ) 
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Observe this system is identical to (7)- (9) with k replaced by* 



sin kAx 

a = — : 

Ax 

Thus the solution to (22) becomes identical to (12), except k 
placed by a , and v by 



(23) 



is re-- 



V = 1 + 



(fe) 



^ sin^kAx 



1/2 



(24) 



Observe that the sinusoidal terms in (12) continue to represent tran- 
sients when V is replaced by v and k by a . Thus the quasi discrete 
case will tend toward a steady state whose transform is: 






Ug(k) = 0 



1 

■ ''o ■"TTTTa ) f "O 



1+X a 



lag y 

— “• h - V 



(25) 



U /I \ U J. H 1 

■ ^o f T ,.2 2 

1+X a 



h - V 
too 



It is easily shown that 



jgg ^ h(x+Ax,0)-h(x-Ax,0) _ 

f o ""o ~^}f 2Ax 



v(x,0) 



(26) 



We shall let 



d(K,t) = v(x,t) _ 



(27) 



Thus if we can Invert 



ia 



2 

1+X a 



and 



.2 2 ’ 

1+X o 
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the steady state solution would be available by convolution. These are 
inverted in Appendix 2, where we show 



- f — 

) l+X 



2 2 
a 



ikx j Ax e 
e dx 



-S|x|/A ^ 

^ 5(x-2nAx), 






and 



00 

fc I 



ia ikx „ -Ax e 
e dk = 



1 u .^2 2 

l+X a 



X'/x^+CAx)" n= 



-e|x|/x 

^ 6 (x-(2n-l)Ax) 



where 



s = s±ntrH¥) . 



, bxim I , 

Ax yx 



(28) 



and 0 < 8 ^ 1 with equality only for Ax = 0. 
Thus we can arrive at 



Vg(x) = v(x, 0) - Ax 



00 

/ 



-3 I s I /X oo 






E 



6 (s-2nAx)( d(x-s,0)ds 



-00 f \X +(Ax) -oo 



1 ^ I 3|2nAx|/X^ 

= v(x,0) ■■■■— \ (e d(x-2nAx, 0) (2Ax) , (29) 

2VX^+(Ax)^ n=-“ 



and similarly 



h (x) = h(x,0) + 

^ / 2 , ,2 
Wx +'■ ' 



IT \ \ ( -S I (2n-l)Ax| /X^ 

" ■ ■■ ■ y Je d(x-(2n-l)Ax,0) (2Ax) 

2XVX‘‘+(Ax)^f n=-“ ( 

(30) 
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Clearly (29)- (30) tend toward the corresponding integral forms (14) as 
Ax 0. 

These expressions can now be examined versus (14) . The conclusions 
that we can draw are that conversion of (l)-(3) to centered, second order, 
quasi-discrete form results in: 

(1) The measure of imbalance, the initial ageos trophic wind, is 
converted to a finite differenced measure. 

(2) Although at steady-state the effects of an initial disturbance 
die off exponentially away from its original neighborhood, the rate of 
decay is less than in the differential case. Furthermore, the rate of 
decay decreases as (Ax/X) increases. 

(3) For the steady states, only the values of the ageostrophic wind 
at alternating points are considered. 

One interpretation of why the sums in (29)-(30) involves only alter- 
nating points may be that, when elimination is tried on (21), one ends 
with the equation for u involving only alternating points: 



at' 



u(x,t) * -f^u(x,t) H ry |u(xf2Ax, t)-2u(x, t)+u(x-2Ax, t)l . 



4 (Ax) 



(31) 



A complete stationary phase analysis of the transient solution in the 
quasi-discrete case: 



u (k,t) = u cos vt + - <v ^ h > 

T o V I o t ° I 



sin vt 



f\, f ^ ^ f ** ) ^ io2 ^ 

v_ (k, t) = — u sin vt + -rr < v ^ h > cos vt 

T M o ,.2)0 f o 



>( 32 ) 



h^(k,t) = 



ioH 'v . _ iaHf j 'v _ ikg v 



u sin vt 
o 



\2 l^o 



£ h / cos vt 
f o 
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is algebraically extremely complicated, although quite straightforward. 
However, the salient features are relatively easily treated, and yield the 
most significant results on the transient behavior. Therefore we only 
present an outline of the details. 

In computing the inversion integrals for (32), terms of the form 

/ 

— CO 

arise. The stationary phase points of these integrals arise as the 
solutions of 



(j)’ (k) = X ± 



^ ^ 2. , , sin kAx 

r A t cos kAx ;; 

Ax 



fl + 



(U 



^ 0 . 



2 2 
sin kAx 



(33) 



This expression can be simplified by adding (-x) to both sides, squaring 

2 2 2 
and writing cos kAx as (1-sin kAx) , to yield a quadratic in sin kAx. 

The quadratic formula then yields as the points of stationary phase the 

solutions of 



sin kAx = 






2.2 ,^2 2 2 .. .2 

A f t -X ) ~4f t X (Ax) 



2 2 2 
2A^f t 



(34) 



It is then easily shown that as t-^ , these solutions closely approximate 



sin kAx = 1 - 



2 2 2 
Aft 



(35a) 



and 



2 2-2 
. 2_ • x^(Ax)^A ^ 

= ■ ,,2,2,2 2, 

(X f t -X ) 



(35b) 
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clearly, for large t , (35b) yields stationary points near both 




and 




(All other solutions of (35b) are beyond the Nyquist cut-off,) The first 
points are slight variants of the stationary points for the differential 
case. The points near ±(7t/Ax) arise solely from the discretization, not 
the physics of the problem. However, since these points give a behavior 
of ((sin kAx)/Ax) identical to that of the stationary points near k = 0, 
it is easily shoxm that they contribute computational transient? with 
precisely the same asymptotic behavior as the physical transients. 

Note, before we consider the effect of terms introduced by the first 
solution, (35a), that a necessary and sufficient condition for the 
stationary phase points to be on the real axis is that the quantity under 
the radical sign in (34) be positive. After some manipulation, this 
condition reduces to 



It is not coincidental that the quantity on the right hand side of the 
inequality is precisely the group propagation velocity for this quasi- 
discrete case. 

Referring again to the contribution from the stationary phase points 
satisfying (35a), observe that we can easily show from (33) that 





t 



> 



t > 0. 
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rck ) = 



//u “ 2 . . 1372 j-(l+(|;)^si»^k4x)sin\Ax + cos\ax 

P(taj ' 



and so, near these stationary points 






±X ft Ax 






Furthermore, note that near these same points 






/(Ax)^+X^ 



(Ax)‘ 
f Ax 



±1 

and a - — 

Ax 



Thus following again the argument of (3.7.5) in Miles [5], we see that the 
points of stationary phase which arise from the solution to the first lead 
to transients whose amplitudes, asymptotically, go as follows: 



r\j A 

U COS vt 

o 



/(Ax) 

2tt X^ft Ax 



1/2 






I ^ 



- h > sin vt“> 
f o ( 



Ax 



2TTX^ft/(Ax)^+X^. 



1/2IX. 

d 



(■ ± ")| 



f 

— u sin vt 



Ax 



,27rX^ft/(Ax)^+ ^ 



1 / 2 , 



“(* to ’ “)| 



<1 V - ft i cos vt^ 

,, 2)0 f o ( 



(Ax)- 



1 1/2IX. 

d 



laH 

— 7 — u sin vt -> 
o 



V 



L2irX^ft 






_2iTX^tAxVAAx) ^+X^ 
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(Note the evaluations are at ±(7r/2Ax) since, as t“X» , these are the 

only solutions of (35a) that also satisfy the Nyquist limit*) 

Viewing the above, it is now clear that the additional stationary 

points which arise in this quasi-discrete case, and which tend toward 

±(tt/2Ax) as t-x» cause two noticeable effects on the transient behavior: 

- 1/2 

(1) All of these transients now die off as t due to dispersion. 

Comparing this to the results in the differential case, we see that these 
transients are at least as persistent as the differential transients, and 

for h(x,t) more so, since the differential transients in h(x, t) die 

^ -3/2 

out as t 

(2) Two of these transients (the first and fifth) are somewhat ill 
behaved as (Ax) 0 . Assuming u(x,0) has only finite power, then the 
Ax term in the denominator should be controlled by tail-off of 

% 7T 

u(± ~ 2 ^ * ® ^ however, these terms are virtually certain to 

be the slowest decaying for small Ax . 

(We note that (35b) also causes stationary points to arise, that tend 
toward ±(tt/Ax) as t-x» , however, the transients from these points do 
not have an asymptotic dispersion decay that depends on (Ax), and decay 
at the same rate as the differential transients.) Although we shall not 
show it analytically, we suspect the additional stationary points arise in 
the quasi-discrete case from two causes, the **folding" in temporal 
frequency that occurs at k = ±(tt/2Ax) , and the high frequency cut off 
at k = ±(tt/Ax) . 
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4 . 



Conclusion 



In this paper we have examined the effect of second order centered 
spatial discretization on a dispersive wave equation. We have shown that 
the methods of Fourier Transforms, and in particular the method of station- 
ary phase, are quite useful in such investigations. In both the differen- 
tial and quasi-discrete cases we have provided closed form expressions for 
the steady-state solutions. These expressions show that the contribution 
from any point in the initial state to the final state decays exponentially 
with distance from that point in both cases, however the rate of decay is 
decreased by discretization. The transients in both cases have been 
analyzed by the method of stationary phase. This analysis shows that the 
discretization introduces stationary phase points that have no counterpart 
in the differential case, and, furthermore, these points contribute 
transients that decay no faster than the differential ones, and in one 
instance, the discrete transient will dominate the differential transient. 
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Appendix 1 



Consider the asymptotic behavior of 



00 



ikx 



A(k)cos vt e dk = \|;(x,t) , 



( 1 " 1 ) 



where V = fvO+x\^ . This integral can be decomposed into 

{ 00 oo 

^ A(k)e^'^l^^^dk + j A(k)e^ 

—OO —00 



♦2<Wdk>, 



where 



(J)j^(k) = kx + vt, and ^2^^^ = kx - vt . 



We can now determine the asymptotic behavior as for fixed x 

of this terra by using the method of stationary phase. Let k^ and k2 
be defined by 

<^[(^ 1 ) = ^2^k2> = 0 . 

But 

2 

- X f 

and so we find as the point of stationary phase for the first Integral: 



k, ^ . (1-2) 

AW:? 



which Is on the path of Integration, for t > x/Xf . Similarly, 
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= 0 yields as the stationary point for the second integral. 



k = 

AW-J 

which is also on the path of integration for t > x/Xf. 
Observe that 



(1-3) 



= -<l>2(k2) = X /Vt^^ , 



and that 



..2^2 2 2,3/2 

<f>"(k.) = -<|.''(kj = = ^± 1 t , ^ 2 > 0 



(Note |f^(k^)| = |<j) 2 (k 2 )| = > for Aft > |x|.) 

Thus we have, by (3.7.5) on p. 51 of Miles [5]> 



-1[ Y + 7 ll 

+ A(k,) e ‘ M 



. f 1 /2^2 2 2 ^ TT 1 

i[ T- f t -X + 7- ] 



11 9 if 2 2 ( i[ Y X^f^t^-x^ + V 

1 P I /A(“k ) e ^ 

2 I r,2^2^2 2.3/2 I ] ^ 2^ 



[A f t -X ]' 



.. 1 J 2.22 2 ^ 7T 

-if T- /X f t -X + T 



+ A( 1 c 2 ) e 



But note that if A(k) is the transform of a purely real valued 
function, 
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A(-k) = 



00 

j a(x)e^^*dk = A (k) . 

“CO 



Thus, we can write: 



'\j 



2 2 
Xf 



il/2 



- ,,2.2^2 2,3/2 

2ir(X f t -X ) 



|r 6 A(k2)cos[ y ^ ] 



+ Im 



* /I N . r 1 /2,2^2 2 . TT J 

A(k2)sin[ Y /X f t -X + "^ ] / » 



(1-4) 



where k 2 is given by (1-3). Observe for Xft > > |xl , the lead term 



acts as 



I 1/2 



.2TrX^ft 



Similarly we can show 



00 

i I 



A(k)sin vt e^^^dk 



% 



■v^2 2 

Xf t 






1/2 



|Re A(k2)sin[ y X^f^t^-x^ + ^ ] 



- Im 



A /I N r 1 /2.2^2 2 . IT , ( 

A(k2)cos[ Y *'1 1 t ^ 1 * 



(1-5) 



Lastly, since the trigonetric functions in both expressions have the 
same frequency, and since 



Xr 6 A(k)]^ + [Im A(k)]^ = |A(k)J 
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we see 



00 


r 2 2 3 


1 A(k)cos vt 'v 

JXoo 


Xf t 


_ ,,2^2^2 2.3/2 

2ir(X f t -X ) 



^ 1/2 



|A(k2)|cos[y 



IT ^ -1, 



- tan [Im A(k. 2 )/Re ACk^)] 

= -| - argCACk^)) 



and 



00 

^ f 

2 tt I 



A(k)sin vt e^^^dx 'v 



,.2 2 
Xf t 



-. 1/2 



2.(x2f2t2-x2)3/2 



|A(k2)|sin[j v^f^t^-x^ + 
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Appendix 2 



Consider 



i. C dx 

2* J 1 + xV’ 



where a 



sin k Ax 
Ax 



(2-1) 



The denominator. 



= 1 + 



(i) 



2 2 
sin k Ax 



( 2 - 2 ) 



is an entire function of k in the complex plane. The zeros of ij^(k) are 
solutions of: 



sin kAx = + 




> 



or, letting k = k^ + i k^ 



sin(k^Ax + i k^Ax) 

= sin (k^Ax) cosh (k^Ax)+i cos (k^Ax)sinh(k^Ax) 




Thus 



sin(k^Ax) cosh(k^Ax) = 0 
cos(k^Ax) sinh(k^Ax) = 



Thus 



k = ± ^ , 

r Ax 



'‘i ■ * fc 



(2-3) 



as opposed to the poles at k^ = 0, k^ = ± 1/X for the continuous case. 
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Thus all zeros of the denominator lie off the real axis. In fact, they are 
specifically distributed as shown: 




Observe further that at its zeros, 



ip* (k) = 2 — sin kAx cos kAx 
Ax 






= ±2l/x^+(Ax)^ 



(2-4) 



where the positive sign holds for k^ > 0, and the negative one for k^ <0. 
Thus the zeros of ij^(k) are simple poles of the integrand. 

Consider (2-1) for x > 0. (Observe that (2-1) is an 
even function of x.) To insure Re(ikx) <0 we will close the contour 
in Im(k) > 0 half plane. As long as we route the contour to avoid all 
the poles there, the integrand is exponentially decaying as |k| , 

Im(k) > 0. Thus, we invoke Jordan’s lemma to yield, for x > 0 



i_ f 

2’ J 1 -f x‘ 



dx 



,2 2 

X a 



= 1 



E 

lm(k)>0 



Residues 



(:*) 



.. E 



ikx 



lm(k)>0 4'’(k) 



k = k + ik. 
r i 
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00 



where 






(Ax)- 



'Yu ‘ to “‘"■’’Hx*)) 



gx 

X 






(Ax)- 




(2-5) 




(2-6) 



Observe (2-5) is not a convergent series in the usual sense* However, 
viewed as a generalized function, it is a Fourier Series for the function, 
periodic of period 2Ax , given in the interval -Ax < x < Ax by 



S(x) == 




where c = 1. 
n 



But observe c 

n 



1 

2Ax 



! 



-Ax 



S(r)e' dr 



1, all n. 



Thus 



S(t) = 2Ax6(t), 



and so its periodic extension becomes 



S(x) = 2Ax 



^ 6(x-2nAx). 

—00 



Thus 
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2it 



CO 




ikx 
e 



2 2 
1+X a 



dk 



Ax 

Vx^+(Ax)^ 



-B|x| 

X 



CO 

^ 6(x-2nAx) . 

—CO 



(2-7) 



Note, if we view 6 as a function 



6(x) 



sinh X , 

X 



it is easily shown: 



6(0) = 1 and 6’ (x) < 0 . 



Thus 0 < 6(x) ^ 1 with equality only at x = 0. 

Since, for Imk > 0 and large 



la 



Ax 



Ax 



T , , 2 2 -V 2 . ^ 2 X Im k 

1+A a X sin kx X e 



the above argument can be essentially repeated for 



CO 

h I 



iae 



ikx 



2 2 
1+X a 



dk = i 



E 



I lg(k)e 



ikx 






-1 



-Bx 

X 



2xVx^+(Ax)^ 



00 

E «>" 



nirx 

Ax 



X > 0 



-1 



2X'^X^+(Ax)‘ 



- g|x| “ • niT / . . ^ 

\ X — > 1 — (x+Ax) 

“ E ■■ 



Ax' 



-(Ax) 



- Bixl 



X\/x\(Ax)^ 



6(x-(2n-l)Ax) 



(2-8) 
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